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Abstract In this paper we investigate a master-slave synchronization scheme of two 
n-dimensional non-autonomous chaotic systems coupled by sinusoidal state error 
feedback control, where parameter mismatch exists between the external harmonic 
excitation of master system and that of slave one. A concept of synchronization with 
error bound is introduced due to parameter mismatch, and then the bounds of 
synchronization error are estimated analytically. Some synchronization criteria are 
firstly obtained in the form of matrix inequalities by the Lyapunov direct method, and 
then simplified into some algebraic inequalities by the Gerschgorin disc theorem. The 
relationship between the estimated synchronization error bound and system 
parameters reveals that the synchronization error can be controlled as small as 
possible by increasing the coupling strength or decreasing the magnitude of mismatch. 
A three-dimensional gyrostat system is chosen as an example to verify the 
effectiveness of these criteria, and the estimated synchronization error bounds are 
compared with the numerical error bounds. Both the theoretical and numerical results 
show that the present sinusoidal state error feedback control is effective for the 
synchronization. Numerical examples verify that the present control is robust against 
amplitude or phase mismatch. 
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1. Introduction 

A lot of efforts have been devoted to the study of synchronization due to its 
wide-scope potential applications in secure communication, 1 " 4 chemical and 
biomedical sciences, " life science, electromechanical engineering " and other 
subjects. Since Pecora and Carroll applied replacing variable control to synchronize 
two chaotic systems, 10 a wide variety of control techniques have been proposed for the 
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synchronization of chaotic systems, such as linear feedback control, " adaptive 
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control, " active control, " sinusoidal feedback control, " and others. " It is 
well known that even though the identical master and slave systems are placed in a 
synchronization scheme, parameter mismatch between two systems often occurs 
because of inevitable perturbation in operations, which can destroy the 
synchronization. " Thus, it is particularly important to investigate robust 
synchronization in the presence of parameter mismatch. 29 " 44 Recently, synchronization 
of a coupled chaotic semiconductor lasers in the presence of parameter mismatch have 
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been studied. " Suykens et al applied standard plant form of synchronization 
scheme to study robust synchronization of master and slave Lur'e systems with 
parameter mismatch via linear state error feedback control, where the master system 
is subjected to an external input. Synchronization of Lur'e systems with parameter 
mismatch was also discussed. " Two-dimensional lattices of diffusively coupled 
chaotic oscillators with parameter mismatch were studied by using the scalar diffusive 
coupling and it was shown that in lattices of Lorenz and Rossler systems the cluster 
synchronization regimes are stable and robust against up to 10%-15% parameter 
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mismatch. It was found in Ref.36 that only negative conditional or transverse 
Lyapunov exponents can not guarantee synchronization of two chaotic systems with 
parameter mismatch. Refs. 37-39 investigated synchronization of chaotic Chua's 
systems in the presence of parameter mismatch. Clearly the contributions mentioned 
above are significant theoretically and practically. However, most of these papers 
concerned synchronization of autonomous systems 29 " 43 and the obtained 

11 11 1-7 1Q 

synchronization criteria were in the form of linear matrix inequalities (LMI), " ' 
most of which must be solved by means of computer algebraic systems, such as 
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Matlab. On the contrary, synchronization criteria in the form of algebraic inequalities 
can reveal the relationship between the criteria and system parameters, so it is more 
convenient to design and analyse the synchronization controller. Synchronization 
criteria for two-dimensional horizontal platform systems with phase mismatch via 
linear state error feedback control was derived in the form of algebraic inequalities. 44 
The authors investigated synchronization criteria of two second-order 
non-autonomous chaotic oscillators coupled by a linear state error feedback control, 
where parameter mismatches may occur in the external excitations or parametric 
excitations. 45 To the best of our knowledge, research on chaos synchronization of 
non-autonomous chaotic systems of higher dimension with parameter mismatch has 
been received little attention up to now. Nevertheless, more and more 
non-autonomous chaotic systems have been found in engineering and physics. 46 " 48 
Especially, high-dimensional non-autonomous chaotic systems are often met with in 
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practice, for examples, a three-dimensional gyrostat system, and a four-dimensional 
loudspeaker system. 46 

In this paper, we investigate a master-slave synchronization scheme of two 
n-dimensional non-autonomous chaotic systems coupled by a sinusoidal state error 
feedback controller, where parameter mismatch exists between the external harmonic 
excitation of master system and that of slave one. One of the advantages of sinusoidal 
controller is that it always produces smooth and bounded output regardless of the 
magnitude of state error variables. Parameter mismatch in this paper implies that 
amplitude, frequency and phase can be all or partly different between the external 
harmonic excitation of master system and that of slave one. 

In the following, we firstly present a master-slave synchronization scheme and 
estimate the synchronization error bound by the Lyapunov direct method and an 
inequality technique. Secondly, the Gerschgorin disc theorem is used to derive some 
synchronization criteria of algebraic inequalities for the scheme. We also establish the 
relationships among the estimated synchronization error bound, parameter mismatch 
and coupling coefficients. The relationships reveal that the synchronization error can 
be controlled as small as possible by increasing the coupling strength or decreasing 
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the magnitude of mismatch. Finally, simulations on a three-dimensional gyrostat 
system verify the effectiveness of the obtained synchronization criteria. The 
difference between the estimated synchronization error bounds and the numerical 
error bounds is also illustrated. Both theoretical and numerical results show that the 
present sinusoidal state error feedback control is effective for the synchronization. 
Numerical examples verify that the present control is robust against amplitude or 
phase mismatch. 

2. Synchronization scheme 

Consider a master-slave synchronization scheme for two non-autonomous 
chaotic systems with different external harmonic excitations coupled by a sinusoidal 
state error feedback controller as follows: 

master system: x = A(t)x + F(x) + m(t) , ( 1 ) 

slave system: y-A(t)y + F(y) + m(t)+u (t) , (2) 

controller: u(t) = (k x sm(x l -y { ), k 2 sin(x 2 - y 2 ),•••, k n sin(x n -y„)) T , (3) 
where the dot represents the derivative with respect to time t , 
x = (x l ,x 2 ,---,x n ) T e R" , y = (y l ,y 2 ,---,y n ) T e R" , A(t) e R" x " is a bounded 

matrix, F is a continuous nonlinear function, m{t) = {m x (t), m 2 (t), ■•■,m n (t)) T e R" is 
external harmonic excitation, m t {i) = f, sin(&>/ + <^ / ) (or m l (t) = f l cos(&>^ + ^ ; ) ), 
l-\,2,---,n . m(t) = (m 1 (t), m 2 (t), ■•■,m n (t)) T eR n , m l {t)=f l sm{d) l t + q) l ) (or 
m,(t) = f l cos(a^t + p l )), f^dtffipj can be different from f^co,,^,, respectively, and 
k l ,k 2 ,...,k n are coupling coefficients. In this paper, we take m,{t) = f t sin(&>^ + <p,) 

and m^t) = f t sin(^ + ^) in the following discussion. The present method can also 

be applied to study the other situations. 

Usually, the synchronization scheme (l)-(3) is difficult to achieve complete 
synchronization due to parameter mismatch. Therefore a concept of synchronization 
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with error bound is introduced as follows. 

Definition 1 The synchronization scheme (l)-(3) achieves synchronization with error 
bound h in a region D if for any bounded initial states of systems (1) and (2) in 
the region/), there exist a constant h >0 and a T > such that _y(f)|</i 

for all t > T , where || • || denote the Euclidean norm. 

According to this definition, our control objective is to select suitable coupling 
coefficients k l ,k 2 ,---,k n such that the distance between trajectories x(t) and y(t) 
of master and slave systems will eventually become less than an error bound h if the 
initial states x(0) and y(0) lie within a certain region D . To this end we need some 
hypotheses. 

Hypothesis 1 There exists a bounded matrix M(x, y) such that 

F(x)-F(y)=M(x,y)(x-y), (4) 

where the elements of M(x, y) are dependent on x and y . 

Fortunately, such a matrix M(x, y) exists in many practical chaotic systems, such 
as Lorenz system, Duffing equation as well as the example in Section 6. 

3 Estimation of synchronization error bound 

Defining an error variable e = x — y , one can obtain the error dynamical system: 

e = A(t)(x -y) + F(x) - F(y) -u(t) + Am = (A(t) - K(t) + M(x, y))e + Am , (5) 
where 

K(t) = diag{k l u l (t), k 2 u 2 (t), • • • , k n u n (t)} , u, (t) = — — , i = 1,2, 

Xi-y t 

Am = m(t) — m(t) = {m { (t) — m y (t), m 2 (t) — m 2 (t), ■ ■ ■ , m n (t) — m n (t)) T e R" , 

m,(t)-m l (t)= f,sm(ca l t + ^) l )-f l sm(a) l t + ^) l ), l-\,2,---,n. 

In what follows Lyapunov direct method is used to analysis the stability of the error 
dynamical system (5) and derive the estimated synchronization error bound. Choose a 
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quadratic Lyapunov function V{e) = e T Se with the positive definite matrix 
S = diag{s l ,s 2 ,...,s n } . Then the derivative of V(e) with respect to time along the 
trajectory of system (5) is 

V{e) = e T [(A(t) - K(t) + M(x, y)) T S + S(A(t) - K{t) + M(x, y))] e + 2Am T Se. (6) 
Obviously, 



ll A HI = \ill(fi sin (^/^ + Pi)-fi sin (*V + ^/)) 2 
/=i 



n n n f + f 

< X \fi sin(^ + tp,)-f l sm(&,t + cp l )\< £(/; +f,) = ^f, • 
1=1 i=i i=i Ji 

Letting G, = ^'~t^' , we have 
fi 

IH|<X/^. 

i=i 

Hence 

2Am T Se<2\\Am T Se\\< 2 ma jAm |||e||< ±f,G,A J|p|, 

/=i 

where A max = max{ 5, . : / = 1,2,- ■■n} is the maximal eigenvalue of matrix 5*. Then 
system (6) can be rewritten as 

V = e T [(A(t)-K(t) + M(x,y)) T S + S(A(t)-K(t) + M (x, y)) + 2 £ />/, ]e - 2 J jjae'e + 2 Am T Se, 

i=i i=\ 

where aeR + . 
Note that if 

|e|>^ with G = max{G,,G 2 , •••,£„}, (7) 
a 

then 

- 2^ 7,ae r e + 2Am r Se < 2\\e\\£ f t (G^ max - a||e||) < . 
i=i i=\ 

Furthermore, we assume that the coupling matrix K(t) can be chosen such that 



Q{t)= 4t(-)Kt iW xy( T £^$ At -(K(±Mt$ + £/ /C ^<2, (8) 

1=1 

which will be discussed in Section 4 in detail. Under the assumptions (7) and (8), we 

know that V <0. Suppose that the balls E(a ) = {e:e <?<a }witha = f 22 - and 

a 

E(a 2 ) = {e : e e < a 2 } are respectively the smallest upper bound and the infimum of the 
ellipsoid E{a { ) = {e : e Se <a { ) . Then we have 

E(a 2 ) = {e: e e < a 2 } zd E(a l ) = {e: e T Se < a { } zd E(a ) = {e : e T e < a } . 
The relationships among E(a ) , E(a 1 ) and E(a 2 ) for a two-dimensional system are 

sketched in Fig.l. V <0 implies that the trajectory of the error system (5) will 
uniformly enter the ellipsoid E(a 1 ) for any trajectory outside the ellipsoid. Hence the 
synchronization scheme (l)-(3) achieves synchronization with an error bound, 
denoted by h , if inequality (8) holds. The radius ^Ja^ of ball E(a 2 ) is an enlarged 

synchronization error bound, namely, h < . Thus H = ^Ja^ is called the 
estimated synchronization error bound, which can be found out in the following. 
Note that S = diag {s l ,s 2 ,...,sj , then 

e Se = s l e l +s 2 e 2 +--- + s a e n . 

Let Z max and /l min be the maximal and minimal eigenvalue of matrix S respectively. 

Hence, the maximal semiaxis i? max and minimal semiaxis i? min of the ellipsoid E{a l ) 
are 



and then 




From the above expression we know that H reaches its minimum if 
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iSj s 2 



■s n -s . In this situation, X n 



I, 



s and the estimated 



synchronization error bound H can be rewritten as 

s 



H = Ga, a = 



a 



(9) 



Now we discuss some special cases of expression (9). If there exist only 
amplitude mismatch and phase mismatch between the external harmonic excitation of 

master system and that of slave one, i.e., a>, = a>, , / = 1, 2, n , then 
f sinO/ + (p,)-f, sm(a>,t + ft ) 

= Wi -ij)sin(e^ + ft) +f l [sm{a) l t + %)- sin(<s^ + ft)]. 
By the differential mean-value theorem, we have 
fsin^t + ft)-^ sinO/ + ft ) 



^ //-// +// \<Pi-Pi\ = fi( 



fi-fi 



\<Pi-&\)- 



As a result, we can replace G, in inequality (7) with 



fx- fx 



+ m-(Pi\ 



Thus 



G = max{G 1 ,G 2 ,--,G„} 



, f fl I ~ I fl fl I „ I fn fn , _ K 

max{ — = — + , ~ — + ~ — + \ ( Pn- ( P„\} 



f 

Especially, if (p l =(p t , then 



f 



fn 



G = max{ 



fl fl fl fl fn fn 



(10) 



(11) 



fl fl fn 

and if f=f,, then 

G = max^ -ft\,\ft -ft\,---,\ft, -ft\}. (12) 

From Eqs.(9)-(12) we know that the synchronization error bound will be sufficiently 
small if the parameter mismatch turns to zero. 



4. Algebraic synchronization criteria 

Now we will deduce the condition under which Eq.(8) holds, that is, the matrix 

Q{t) is negative definite. Assuming that 

{A{t) + M{x,y)) T S + S{A(t)+M{x,y)) = {sb lj {t,x,y)) nxn , 
where b ij {t,x,y)-b ji {t,x,y) , i,j-\,2,...,n, then the matrix Q{t) mentioned in 
expression (8) becomes 

sb n {t,x,y)-2k l su l +lY J fia sb n (t,x,y) ... sb ln (t,x,y) 

1=1 

sb 2l (t,x,y) sb 22 (t, x, y) - 2k 2 su 2 + 2^ f t a ... sb 2ll (t,x,y) 



2(0 = 



sb nl (t,x,y) sb n2 (t,x,y) ... sb m (t,x,y) -2k n su n +2^f l a 



V 



Our object is to choose coupling coefficients k t such that the matrix Q(t) is 
negative definite. To this end, we first introduce the Gerschgorin disc theorem. 
Lemma 1. (Section 6.1 in Ref.49) LetZ = (z..)ei?" x "and p l ,p 2 ,---,p„be positive 
numbers. Then all the eigenvalues of Z lie in the region 

u{A e C:|A-z,|<^£^l4> 



Pi 7 =1 

J*' 



where C is the set of complex numbers. 

According to Lemma 1, the matrix Q(t) is negative definite if k^p^a and s 
satisfy 



2k i su i > sb u (t, x, y) + 2^ f t a + — ^ Pj \b tj (t, x, y)\, i = \,2,...,n, (1 3) 

i=i Pi j=\ 

j*i 



namely 



2k i u i >b ii (t,x,y) + g — + — ^Pjfa(t,x,y)\, z=l,2,...,n. (14) 

H Pi 7=1 

j*i 

According to Definition 1, if the slave system is synchronized with the master with 
error bound h ( h<H ), then x i - H < y i < x i + H,i = 1,2, • • • , n . Recall that the 



elements of A(t) and M (x, y) are bounded, there exist constants by > such that 
\by (t, x,y)\< by , i^j, and constants b u such that b H (t,x,y)<b u , i - 1, 2, .. ., n , for a 
prescribed error bound H . Substituting the bounds of by(t,x, y) and b u (t,x, y) into 
the inequalities (14), we obtain more rigorous inequalities as Hypothesis 2 

2G ifi j 

2kn>b u + + — ipjby, i = \,2,...,n. (15) 

H Pi 7=1 

If the coupling coefficients k i and parameters p t are chosen such that inequalities 

(15) hold, then the matrix Q(t) is negative definite, which implies that the trajectory 

of the error system (5) will uniformly enter the ellipsoid E(ci x ) (so the error bound 

h < H) for any trajectory outside of the ellipsoid E(ci x ) . We have therefore outlined a 
proof of the following theorem. 

Theorem 1 The synchronization scheme (l)-(2) achieves synchronization with error 
bound h < H if the Hypotheses 1 and 2 are satisfied. 

In the following, we will derive some easily verified synchronization conditions 

from Theorem 1 . Note that u t (t) > , i - 1,2, • • • , n , if (xj , x 2 , ■ ■ • , x n ) and (>>, , y 2 , • • • , y n ) 

are limited in the region D = —y t \ < n, i = 1,2,- • -,n} . Under this restriction, we 
can rewrite inequalities (15) as 

k i >—(b ii +^— + -Y.PM i = l2,-,n. (16) 

iu i n p i J=l 

■ In 

If the variables are further restricted in the region D = {\x i —y i \<—,i = l,2,---,n}, 

2V2 

then we have 1 ^ » / = 1,2, . As a result, simpler criteria can be 

obtained from inequalities (16) as follows: 
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Zn 2G ^' 1 - 

*' > ^ + ^~ + ^* '- = l>2,..,n. (17) 

We will apply the criterion (17) instead of (15) or (16) in the upcoming example. 



5. Analysis of synchronization error bound 

In this section we will discover some relationships among the estimated 
synchronization error bound, coupling coefficient and parameter mismatch. Assume 

that there exists an i e {1,2, •••,«} such that 

max{^= (b a + + — f J P i b ii )}=^ i + + — 1> A / ) • 

For simplicity of notation, in what follows we denote i as i . If the coupling 

coefficients k i are chosen identically, i.e., K{t) = diag{ku x {t),ku 2 {t),--- ,ku n (t)} , 
then the synchronization criteria (17) become 

In 2G & 1 1 - 

k> TTr ib « + ^— + -^Ph ) > ^{i' 2 '-'"}- (18) 

In practice, we can take 

where £ is a small positive constant. For the case of only amplitude mismatch in the 
external harmonic excitations, substituting expression (11) into the above inequalities 
yields 



371 



2Zfi\f-f 



k --^ + ^T^ + p-p A)+£ - (20) 

Hence we obtain the relationships among the estimated synchronization error bound 
H , the amplitude mismatch / - / and the coupling coefficient k , 
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6*\f-f\Zf' 



H = 



1=1 



Pi m 



(21) 



For a fixed k , F^O as / - / 



->0 



Similarly, for the case of only phase mismatch in the external harmonic 
excitations, substitution of expression (12) into expression (19) yields 



3k A<P~<P^Lfl , n 

k = (b u + + — Yp A ) + e 



(22) 



Hence we obtain the relationships among the estimated synchronization error bound 
H , the phase mismatch \<p - (p\ and the coupling coefficient k , 



H = 



i=i 



4j2(k-e)-3xb H - — f jPj b ij 



(23) 



For a fixed k, H — >0 as \(p-(p\^>-0 . 



6. Numerical example 

Before starting with the example, we first outline the design procedure as follows: 

1) Construct the slave system associated with the master system via sinusoidal state 
error feedback control. 

2) Verify Hypothesis 1 to be satisfied. 

3) For a prescribe error bound//, determine the inequalities (17) instead of 
Hypothesis 2 to ensure that the synchronization scheme (l)-(3) achieves 
synchronization with error bound h<H . 

4) Analyse the relations among the estimated error bound, coupling coefficient and 
parameter mismatch. 

5) Verify the control robustness against particular parameter mismatches. 

22 

Consider a gyrostat system 
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Xj = -x 2 x 3 - 0.5(1 + 6.5 cos t)x 2 + 0Ax 3 - 0.002(xj + xf) 

x 2 =x x x 3 -0.4x 3 +0.5(1 + 6.5 cos -0.002(x 2 +x\) . (24) 
x 3 = -0.2xj + 0.2x 2 - 0.2 x 3 - 0.002(x 3 + x] ) + 1 .625 sin t 

Obviously, m(t) is (0, 0, m 3 (t)) T with m 3 (t) = f 3 sin(&> 3 £ + q> 3 ) = 1 .625 sin t , compared 

with system (1). The gyrostat system exhibits chaos behavior, as shown in Fig. 2(refer 

to Ref.22 for more details). Subjected to a sinusoidal error feedback control, the slave 

system is constructed as follows: 

X = -y 2 y3 - °- 5 (l + 6.5 cos t)y 2 + 0.4j; 3 - 0.002(7! + y\) + k x sin^ - y x ) 

y 2 = y 1 y 3 -0Ay 3 +0.5(l + 6.5cost)y l -0.002(y 2 +y 3 2 ) + k 2 sm(x 2 -y 2 ) . (25) 

y 3 = -0.2y l + 0.2y 2 - 0.2y 3 - 0.002(>> 3 +y 3 3 ) + f 3 sm(cb 3 t + q> 3 ) + k 3 sin(x 3 -y 3 ) 

The error variable e = (x l - y l , x 2 - y 2 , x 3 - y 3 ) T satisfies 



e = (A(t) - K{t) + M{x, y))e + Am 



with 





( -0.002 


-0.5(1 + 6.5 cos t) 


0.4 ' 




( k x u x N 


A(t) = 


0.5(l + 6.5cosf) 


-0.002 


-0.4 


, *(0 = 


k 2 u 2 




v -0.2 


0.2 


-0.202; 




v k 3 u 3j 


u,(t) = 


sin(x,-j,) . 

> i 

*, - y. 


1,2,3, 









Mix,y) = 



-O.OOlixf + x x y x +y 2 x ) 

o 



+. 
o 



-x, 



-0.002ix 2 2 +x 2 y 2 +y 2 2 ) 



-0.002ixl +x 3 y 3 +y 3 ) 



Hence 

Qit) = (A(t) - K{t) +Mix, y)f S + S(A(t) - Kit) +M(x, y)) + 2/ 3 a/ 2 







s(-0.004 - 0.004(x 2 + x x y x + y 2 x )) 
-2sk x u x + 2 f } a 

s(-0.004 - 0.004(x 2 2 + x 2 y 2 + y 2 2 )) 
-2sk 2 u 2 + 2 f 3 a 

s(-0.404-0.004(x 2 +X3J3 +y] )) 
-2s£ 3 m 3 + 2f 3 a j 







s(0.2-x 2 ) 



s(0.2-x 2 ) 
si-0.2 + x x ) 



s(-0.2 + x,) 



Under the assumption that the slave system can be synchronized with the master one 
with error bound h ih <H), there exists a T > such that \x t —y i \ <H,i = 1,2,3 , 
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for t > T . Hence, 

b n (t,x,y) = -0.004 -0.004(x 2 +x x y l + y 2 ) 

= -0.004x, 2 - 0.004^ y, - 0.004 j, 2 - 0.004 
<-0.004jc 1 (x 1 ±i/)- 0.004 
= -0.004(x 1 ± y) 2 + 0.001// 2 - 0.004 
< 0.00 IH 2 -0.004. 

The estimated synchronization error bound H can be chosen according to the 
demand in practical application. In this paper, it is set to be < H < 1 , then we have 

b u (t, x, y) = -0.004 - 0.004(x 2 + x x y x +yf)< 0.001 - 0.004 = -0.003 . 
Similarly, 

b 22 (t, x, y) = -0.004 - 0.004(x 2 + x 2 y 2 + y 2 2 ) < 0.00 1 - 0.004 = -0.003 , 

b 33 (t, x, y) = -0.404 - 0.004(x 2 + x 3 y 3 + y 2 ) < 0.00 1 - 0.404 = -0.403 . 

Therefore, the synchronization criteria corresponding to inequalities (17) can be 
obtained as follows: 

k >-^= (-0.003 + ^G+^ 10.2 -x 2 1), 
4V2 H pi 

k 2 >^ (-0.003 + ^-G + ^1-0.2 + x, I), (26) 
4V2 H p 2 

k 3 > ^= (-0.403 + ^ G + 2l | .2 - x 2 1 + ^ 1-0.2 + x, I) . 
4V2 # P 3 ^3 

Thus we conclude that the synchronization scheme (24)-(25) can achieve 

synchronization with error bound h < H if the coupling coefficients k x , k 2 and k 3 

satisfy the inequalities (26). 

From Fig.2, we know that the bounds of the chaotic attractor are -3 < x x < 2.4, 

-3 < x 2 < 2 and-2 < x 3 < 1.4 . If we take f 3 = 1.6 , a> 3 = 0.9 and q> 3 = 0.1 , namely, 
the amplitudes, frequencies and phases of external harmonic excitations of two 

14 



coupled systems are mismatched simultaneously, then G ■ 



A+A 



2.02. For a 



A 

prescribed estimated synchronization error bound, such as H = 0. 1 , and p l = p 2 = p 3 , 
we can get ^> 113.288, k 2 > 1 13.288 and k 3 > 117.964 from criteria (26). The 

simulation in Fig. 3 shows that the synchronization scheme (24)-(25) achieves 
synchronization with error bound h < H = 0.1, where the coupling coefficients are 

^ = 1 13.29 , k 2 =113.29 and k 3 = 117.97 , and the initial values are 
(x, (0), x 2 (0), x 3 (0)) = (1,1,1) and (y, (0), y 2 (0), y 3 (0)) = (-1, - 1, - 1) . 

If the frequencies of two coupled systems are same, we take / 3 =1.6, 

f.-f. 

cl) 3 = co 3 = land q) 3 = 0.1 , then we can get G = — = \-\<p 3 -(f> 3 \ = 0.12 . If we choose 

A 

H = 0.01 and p l =p 2 =p 3 , we can get k x > 69.467 , k 2 > 69.467 and k 3 > 74.143 

from criteria (26). As shown in Fig. 4, the synchronization scheme (24)-(25) can 
achieve synchronization with error bound h<H = 0.01, where the coupling 

coefficients are k x = 69 .47 , k 2 = 69 .47 and k 3 = 74.15 , and the initial conditions are 

the same as in Fig. 3. 

If there exists only amplitude mismatch between two coupled systems, we take 

A- A 



f 3 = 1 .6 , a> 3 = a> 3 = 1 and <p 3 = <p 3 = , then G 



0.02. If H = 0.0l and 



A 

Pi = p 2 = p 3 > we can get k x > 16.027 , k 2 > 16.027 and k 3 > 20.703 from criteria 
(26). Fig. 5 shows that the synchronization scheme (24)-(25) achieves synchronization 
with error bound h<H = 0.01, where the coupling coefficients are ^=16.03, 

k 2 = 16.03 and k 3 = 20.71 , and the initial conditions are the same as in Fig. 3. 

/ -f 

In the case of & 3 = co 3 and q> 3 =(p 3 , we know that G = — = from expression 

A 

(11). Taking p x = p 2 = 100 and p 3 =135 , another criteria for synchronizing systems (24) 
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and (25) can be obtained from inequalities (26) 

A- f 3 



In 2 
h >^= (-0.003 + 

4^2 



H 

fs-fs 



3k 2 
k 2 >^= (-0.003 + - 

4^2 H 

7 3^ , _ 2 /j - / 3 

fc 3 > — ^= (-0.403 + ■ 



135, 
100 1 1 



1 35 

+ — 1-0.2 + xJ), 
100 1 1 



100 lno , 100, no ,, 

+ 0.2 -x 9 + -0.2 + x ). 

135 1 1 135 1 



(27) 



4>/2 v H 

Substituting the bounds of the strange attractor of the gyrostat system into inequalities 
(27), we obtain more rigorous inequalities as follows 



k x > 7.209 + - 



3.34 



fs-fs 



H 



k 2 > 7.209+ 



3.34 



A- A 



H 



(28) 



L > 7.244 + - 



3.34 



/ 3 -/ 3 



7/ 



The following criterion is corresponding to criteria (20) 



£ = it =7.244 + - 



3.34 



/ 3 -/ 3 



7/ 



+ £ . 



(29) 



For simulation, £ = 0.001 is chosen. The estimated error bound H corresponding 
to expression (21) becomes 



H = 



3.34 



/ 3 "/ 3 



£-7.245 



Taking / 3 -/ 3 = 0.02 as an example, then 



(30) 



7/ = 



0.0668 



(31) 

it- 7.245 

For different coupling coefficients k , the estimated synchronization error 
bounds (31) are compared with numerical error bounds, which is shown in Fig. 6. The 
numerical synchronization error bound is calculated by 

||e|| max = max ^/(xj -y x ) 2 + (x 2 -y 2 ) 2 +(x 3 - y 3 ) 2 , where T Q represents a time 

threshold after which the synchronization error will be stabilized in a prescribed error 
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bound and T is a sufficiently large positive number. In the following simulations, 



T Q = 200 and T = 10000 time unit are taken. 



For physical meaning, we calculate the relationship between the relative error 



bound H(%) = — — and relative parameter mismatch Af (%) = 

w / 

II Umax J 
3.34/3 



as follows: 



H(%) = ^ '—^ Af 3 (%) , (32) 

x (A: -7.245) 3 



where Ibcll is defined as lljcll = max Jx, 2 + x, 2 +x, 2 . For the gyrostat system (24), 

II Umax II Umax T <t<T 

IWLax = • Taking k = 15 as an example, then 

//(%) = O.2A/3 (%). (33) 

Let \\e\\ (%) = jj (( max . The estimated synchronization error bounds H(%) are 

II Umax L/C 

II Umax 

compared with the numerical error bounds ||e|| (%) versus different amplitude 
mismatches A/ 3 (%) . The comparison result is shown in Fig.7, from which we know 

that a 12% amplitude mismatch induces approximately only a 2.4% synchronization 
error, which shows that the present control is robust against amplitude mismatch. 

If there exists only phase mismatch between the two coupled systems, we take 

f 3 =f 3 =1.625. ft> 3 =ft> 3 =1 and^ 3 =0.1 . Hence G = \<p 3 -<p 3 \ =0.1 . We can get 
k x > 59.614, k 2 > 59.614, k 3 > 64.29 from criteria (26) if H = 0.01 and p x = p 2 = p 3 . 
Fig. 8 shows that the synchronization scheme (24)-(25) achieves synchronization with 
error bound h<H = 0.01 . The coupling coefficients are k x = 59.62 , k 2 = 59.62 and 

k 3 = 64.3 , and the initial conditions are the same as in Fig. 3. 

For the case of f 3 = f 3 ,&> 3 = a> 3 , from expression (12) we know that G = \cp 3 —<p 3 \ ■ 
If p l =p 2 =\00, p 3 =\35, another criteria for synchronizing the systems (24) and 
(25) can be obtained from inequalities (26) 



17 



^ > ^ (-0.003 + 2 ^3 1^3 ~ ^ 3 1 + 135 _ 

1 4^f2 H 100 1 21 

k 2 ^(-0.003 + ^^ + ^1-0.2 + xJ), (34) 

2 4^ H 100 1 11 

2>7i IfAm.-mA 100, ,100, 

L > ^= (-0.403 + ^ + — 0.2 - x 2 + — -0.2 + x ) . 

4>/2 # 135 1 21 135 1 11 

Substituting the bounds of the chaotic attractor of the gyrostat system into inequalities 
(34) leads to more rigorous inequalities 

5.428 I0, -0,1 
h > 7.209 + ^ ^ 3| , 

1 H 

5.428 I0, -0,1 

> 7.209 + ^ 3| , (35) 

2 H 

5.4281©, -0,1 

> 7.244 + ^ 3| . 

3 // 

The following criterion is corresponding to criteria (22) 

5.428 10, -0,1 

£ = =7.244 + |r3 r3| +g. (36) 

3 H 

In simulation, £ = 0.001 is chosen as above. The estimated error bound H 
corresponding to expression (23) becomes 



5.428 10, -0,1 

# = ^ ^ 3| . (37) 

£-7.245 

Taking |0 3 - #? 3 1 = 0. 1 as an example , then 

H = ^^- . (38) 
£-7.245 

For different coupling coefficients k , the estimated synchronization error bounds (38) 
are compared with the numerical error bounds, shown in Fig. 9. 

I0, - 0, 1 

Similar to the discussion of Eq.(32), letting A0 3 (%) = - '-, we can obtain the 

3 

relative error bound as follows: 
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H{%) = ,, 5A2S(P3 A^(%) . (39) 

x (A: -7.245) 3 

II Umax v ' 

Taking cp 3 = 2n and k = 50 as an example, then 

//(%) = 0.228Ap 3 (%). (40) 
For different phase mismatches A ^ 3 (%) , the estimated synchronization error bounds 

(40) are compared with the numerical error bounds ||e|| (%), shown in Fig. 10. 

From the comparison we know that a 12% phase mismatch induces approximately a 
2.7% synchronization error. So we conclude that the present control is robust against 
phase mismatch. 

7. Conclusions 

In this paper we have investigated in detail the synchronization of a class of 
master-slave non-autonomous chaotic systems coupled by a sinusoidal state error 
feedback controller, where parameter mismatch exists between the external harmonic 
excitation of master system and that of slave one. A concept of synchronization with 
error bound is introduced. The synchronization criteria are in the form of algebraic 
inequalities, so they are convenient in applications. The bounds of synchronization 
error are estimated analytically. It follows from the estimated synchronization error 
bound that the synchronization error can be controlled as small as possible by 
increasing the coupling strength or decreasing the magnitude of mismatch. The 
three-dimensional chaotic gyrostat system is used as an example to verify the 
effectiveness of the proposed synchronization criteria. The comparisons of estimated 
synchronization error bounds with the numerical error bounds show that a 12% 
amplitude mismatch or phase mismatch induces acceptable small synchronization 
errors, which verifies that the present control is robust against amplitude or phase 
mismatch. 
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Figure captions 

Fig. 1 Illustration of circles E(a ), E(a 2 ) and ellipse E{a l ) 

Fig.2 Phase portrait of the gyrostat system (24) showing the strange attractor 

Fig. 3 Error between the master-slave systems (24)-(25) with =113.29 > k 2 =113.29 

and k 3 =117.97 when amplitude, frequency and phase are mismatched simultaneously. 

The beeline and curve represent the estimated error bound H = 0. 1 and numerical error 
respectively. 

Fig.4 Error between the master-slave systems (24)-(25)with k x = 69 .47 , k 2 = 69 .47 
and ^ 3 = 74.15 when amplitude and phase are mismatched. The beeline and curve 
represent the estimated error bound H = 0.01 and numerical error \e(t)\ , respectively. 
Fig. 5 Error between the master-slave systems (24)-(25) with k x =16.03 , k 2 =16.03 
and k 3 =20.71 when there exists only amplitude mismatch. The beeline and curve 
represent the estimated error bound H = 0.01 and numerical error \e(t)\ , respectively. 

Fig. 6 Synchronization error bounds versus the different coupling coefficients k . The 
solid curve and dotted curve represent the estimated synchronization error bound H 
and the numerical error bound llell , respectively. 

II Umax A * 

Fig. 7 Synchronization error bounds versus the different amplitude mismatches A/ 3 (%). 
The solid curve and dotted curve represent the estimated synchronization error bound 
//(%)and the numerical error bound || e|| m ax( 0// °)' res P ect ively. 

Fig. 8 Error between the master-slave systems (24)-(25) with =59.62, k 2 = 59.62 
and k 3 = 64.3 when there exists only phase mismatch. The beeline and curve 
represent the estimated error bound H = 0.01 and numerical error \e(t)\ , respectively. 

Fig. 9 Synchronization error bounds versus the different coupling coefficients k . The 
solid curve and dotted curve represent the estimated synchronization error bound H 

24 



and the numerical error bound I el , respectively. 

II Umax * J 

Fig. 10 Synchronization error bounds versus the different phase mismatches A<p 3 (%) . 

The solid curve and dotted curve represent the estimated synchronization error bound 
H{%) and the numerical error bound II ell (%), respectively. 
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Fig.l Illustration of circles E(a ) and E(a 2 ) , and ellipse E(a 1 ) 
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Fig. 2 Phase portrait of the gyrostat system (24) showing the chaotic attractor 
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Fig.3 Error between the master-slave systems (24)-(25) with ^ =113.29 » & 2 =113.29 and 
fc 3 =117.97 when amplitude, frequency and phase are mismatched simultaneously. The beeline 
and curve represent the estimated error bound// = 0. 1 and numerical error \\e(t)\\ , respectively. 
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Fig.4 Error between the master-slave systems (24)-(25) with k x = 69 .47 , k 2 = 69 '.47 and 
^ 3 = 74. 15 when amplitude and phase are mismatched. The beeline and curve represent the 
estimated error bound H = 0.01 and numerical error \\e(t)\\ , respectively. 
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Fig.5 Error between the master-slave systems (24)-(25) with ^ =16.03 > k 2 =16.03 and 
k 3 =20.71 when there exists only amplitude mismatch. The beeline and curve represent the 
estimated error bound H = 0.01 and numerical error ||e(0|| > respectively. 
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Fig.6 Synchronization error bounds versus the different coupling coefficients k . The solid curve 
and dotted curve represent the estimated synchronization error bound H and the numerical error 
bound Ml , respectively. 
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Fig.7 Synchronization error bounds versus the different amplitude mismatches A/ 3 (%) . The solid 
curve and dotted curve represent the estimated synchronization error bound H(%) and the 
numerical error bound II ell (%), respectively. 
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Fig.8 Error between the master-slave systems (24)-(25) with = 59.62 , k 2 = 59.62 and 
£ 3 = 64.3 when there exists only phase mismatch. The beeline and curve represent the estimated 



error bound// = 0.0 land numerical error ||e(/)L respectively. 



29 



. Oft 



o 
I 



0.06 



0.04 



0.02 



14 



16 



18 20 
k 



22 



14 



Fig. 9 Synchronization error bounds versus the different coupling coefficients k . The solid curve 
and dotted curve represent the estimated synchronization error bound Hand the numerical error 

bound e , respectively. 
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Fig. 10 Synchronization error bounds versus the different phase mismatches A #> 3 (%) • The solid 
curve and dotted curve represent the estimated synchronization error bound H(%) and the 
numerical error bound e (%)> respectively. 
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